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With  the  recent  advances  in  short  range  infrared  guided  missiles,  a  signi¬ 
ficant  change  in  close  air  combat  characteristics  emphasizing  instantaneous 
maneuverability  becomes  inevitable.  This  has  led  to  a  growing  interest  in  the 
so  called  "supermaneuverability"  of  fighter  aircraft,  Where  previously 
unattained  regions  of  the  maneuver  envelope  are  attempted.  For  example,  post¬ 
stall  flight  maneuver  and  side  slipping  (Herbst,  1981;  1983).  There  exists  a 
need  for  improving  existing,  and  devising  novel,  vortex  control  devices  that 
might  help  achieving  the  desired  high  degree  of  maneuverability. 

In  a  recent  research  program  conducted  at  Flow  Research,  it  was  discovered 
that  the  leading  edge  vortex  on  a  delta  wing  at  a  constant  angle  of  attack 
consists  of  a  series  of  discrete  smaller  vortices.  These  vortices  are  shed  at 
a  predictable  frequency.  Their  trajectories  follow  the  general  outline  of  the 
classical  large  bound  vortex  on  each  side  of  the  wing.  As  the  discrete 
vortices  follow  this  course,  two  vortices  would  begin  to  roll  around  each 
other  and  merge  to  form  a  single  larger  vortex.  This  process  was  observed 
to  repeat  itself  several  times. 

Preliminary  attempts  to  modulate  the  shedding  and  the  pairing  of  the 
discrete  vortices  were  successful.  The  flow  at  the  leading  edge  of  the  delta 
wing  was  perturbed  by  impulsively  injecting,  at  a  given  frequency,  a  secondary 
fluid  from  a  slot  along  that  edge.  The  flow  field  was  significantly  affected 
when  the  injection  frequency  was  chosen  to  be  a  subharmonic  of  the  shedding 
frequency  of  the  discrete  vortices. 

The  ultimate  goal  of  the  present  work  is  to  develop  a  device  that  can  be 
used  to  change  the  shedding  and  the  pairing  of  the  discrete  vortices  on  a 
delta  wing  aircraft.  The  device  will  be  optimized  to  significantly  change  the 
lift  of  the  airplane  at  constant  speed  and  constant  angle  of  attack.  By 
applying  the  method  to  both  leading  edges,  the  total  lift  of  the  wing  will  be 
altered;  alternatively,  by  using  the  perturbation  preferentially  on  only  one 
aide  of  the  wing,  the  rolling  moment  around  the  axis  of  symmetry  of  the 
aircraft  is  controlled.  Such  a  device,  when  successfully  developed,  will 
enable  the  pilot  of  a  delta-wing  fighter  aircraft  to  achieve  a  previously 
unattained  degree  of  maneuverability. 


During  Phase  I  of  this  research  program,  experiments  were  conducted  in  a 
water  towing  tank  and  a  wind  tunnel  to  characterize  the  flow  field  around  a 
delta  wing  with  and  without  the  vortex  control  device.  Flow  visualizations, 
velocity  surveys  using  fast-response  probes  and  aerodynamic  force  measurements 
were  conducted.  In  this  report,  a  brief  background  is  given  in  Section  2. 

The  experimental  approach  is  described  in  Section  3.  The  flow  visualization, 
velocity  and  force  measurements  results  are  detailed  in  Sections  4  through  8. 
Finally,  a  summary  and  recommendation  for  future  research  are  found  in  the 
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The  flow  over  delta  wings  at  constant  angle  of  attack  is  dominated  by  two 
large  bound  vortices  that  result  from  the  flow  separation  at  the  leading  edge 
(Bile,  1958;  Wentz  &  Kohlman,  1968;  Werle,  1973;  Koerner  &  Borst,  1975;  and 
Faery  et  al.,  1981).  With  a  sharp  leading  edge  at  an  angle  of  attack  a,  the 
flow  is  separated  along  the  entire  leading  edge  forming  a  strong  shear  layer. 

The  shear  layer  is  wrapped  up  in  a  spiral  fashion  resulting  in  a  large  bound 
vortex  on  each  side  of  the  wing  as  sketched  in  Figure  la.  The  two  vortices 
appear  on  the  suction  surface  of  the  wing  in  the  form  of  an  expanding  helix 
when  viewed  from  the  apex  of  the  lifting  surface.  The  low  pressure  associated 
with  the  vortices  produces  additional  lift  on  the  wing,  often  called  nonlinear 
or  vortex  lift,  which  is  particularly  important  at  large  angles  of  attack. 

A  recent  research  program,  supported  by  the  Air  Force  Office  of  Scientific 
Research,  was  conducted  at  Flow  Research  to  determine  and  examine  the  structure 
of  the  shear  layer  at  the  leading  edge  of  a  delta  wing  and  other  three- 
dimensional  lifting  surfaces  in  steady  and  unsteady  flights  (Gad-el-Hak  et  al., 
1983a;  1983b;  1985b;  1985c;  1985d;  1986a;  1986b;  1986c).  Flow  visualization 
experiments  clearly  showed  that  the  shear  layer  flow  near  the  leading  edge  of  a 
delta  wing  at  constant  angle  of  attack  rolls  up  into  discrete  vortices  as 
sketched  in  Figure  lb.  The  trajectory  of  these  vortices  followed  the  general 
outline  of  the  large  bound  vortex.  As  they  followed  this  course,  two  vortices 
would  begin  to  roll  around  each  other  and  merge,  or  pair,  to  form  a  single 
layer  vortex,  much  the  same  as  observed  in  a  free  shear  layer  (Figure  2) 
originating  at  a  splitter  plate  between  two  streams  of  differing  velocity 
(Winant  &  Browand,  1974;  Brown  &  Roshko,  1974;  and  Roshko,  1976). 

Examples  of  the  discrete  vortices  on  a  delta  wing  in  steady  flight  are  shown 
in  Figures  3  and  4.  In  Figure  3,  a  60*  sweep,  sharp-leading-edge  delta  wing  was 
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towed  in  an  18-meter  water  channel  at  a  root-chord  Reynolds  number  of  1.3x10  , 

and  at  an  angle  of  attack  of  10*.  Fluorescent  dye  was  injected  from  a  slot 

near  the  left  hand  side  leading  edge  and  illuminated  with  a  bank  of  black 

lights.  The  dye  marks  the  discrete  vortices  with  their  axes  emanating  from  the 

apex  of  the  delta  wing.  Figure  4  is  a  side  view  of  a  45*  sweep  delta  wing 

Which  had  a  BACA  0012  profile  at  each  spanwise  section.  The  wing  had  a 

4 

root-chord  Reynolds  number  of  2.5x10  .  The  fluorescent  dye  is  illuminated 
with  a  vertical  sheet  of  laser  light  (Gad-el-Hak  et  al.,  1981)  parallel  to  the 
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flow  direction  at  z  =  10  cm  (40H  of  the  root-chord  length).  In  this  photograph, 
the  pairing  process  between  two  discrete  vortices  is  clearly  indicated. 

For  angles  of  attack  below  5*.  the  pairing  process  seemed  to  be  inhibited 
by  the  proximity  of  the  vortices  to  the  wing  surface.  The  discrete  vortices 
were  observed  on  delta  wings  with  sharp  and  with  blunt  leading  edges,  with  delta 
wings  having  different  sweep  angles,  and  with  a  Reynolds  number  varying  over  at 
least  one  order  of  magnitude.  The  shedding  frequency  of  the  discrete  vortices 
depends  on  the  mean  velocity  of  the  wing  and  the  parameters  effecting  the  thick¬ 
ness  of  the  shear  layer.  This  frequency  can  be  modulated  by  a  pitching  motion 
of  the  wing  (Gad-el-Hak  et  al.,  1983b). 

It  should  be  noted  that  the  shedding  of  discrete  vortices  in  a  free  shear 
layer  has  been  observed  to  occur  resonantly  at  all  Reynolds  numbers  investigated 
in  different  laboratories  around  the  world  (Ho  &  Huerre,  1984),  including 
Reynolds  numbers  well  in  excess  of  the  critical  one.  We  believe  that  the  analo¬ 
gous  phenomenon  of  shedding  from  a  delta  wing  is  also  occurring  at  all  Reynolds 
numbers  and,  hence,  our  observations  have  universal  validity.  The  shedding  of 
discrete  vortices  from  a  free  shear  layer  or  from  the  leading  edge  of  a  delta 
wing  is  to  be  contrasted  to  the  shedding  from  a  bluff  body.  In  the  latter 
case,  the  shedding  occurs  resonantly  at  subcritical  Reynolds  numbers  but  occurs 
more  randomly  at  higher  Reynolds  numbers. 

During  the  above  mentioned  research  program,  attempts  were  made  to  modulate 
the  shedding  and  the  pairing  of  the  discrete  vortices  from  a  delta  wing  at 
constant  angles  of  attack.  This  was  inspired  by  the  recent  attempts  to  control 
the  growth  rate  of  a  free  shear  layer  by  perturbing  the  flow  at  the  end  of  the 
splitter  plate  (Ho  &  Huerre,  1984).  The  flow  at  the  leading  edge  of  the  delta 
wing  was  perturbed  by  impulsively  injecting  a  secondary  fluid  from  a  slot  along 
that  edge.  The  frequency  of  perturbation  was  varied  over  a  wide  range,  and  it 
was  found  that  maximum  changes  in  the  flow  field  occurred  when  this  frequency 
was  a  subharmonic  of  the  natural  shedding  frequency  of  the  unperturbed  wing. 


3.  gTPKBTHCTTAL  APPROACH 


3.1  Towing  Tank 

The  present  experiments  were  conducted  in  both  a  water  towing  tank  and  a 

high-speed  wind  tunnel.  The  flow  visualization  tests  were  best  suited  for  the 

water  facility.  The  towing  tank,  Figure  5,  is  18  meter  long,  1.2  meter  wide, 

and  0.9  m  deep  (see  Gad-el-Hak  et  al.,  1981).  The  delta  wing  was  rigidly 

mounted  on  a  carriage  that  rides  on  two  tracks  mounted  on  top  of  the  towing 

system.  During  towing,  the  carriage  was  supported  by  an  oil  film.  Which  ensured 

a  vibrationless  tow.  The  equivalent  freestream  turbulence  was  checked  using  a 

hot-film  probe,  and  was  found  to  be  about  0.1H  of  the  mean  towing  speed.  The 

carriage  was  towed  by  two  cables  driven  through  a  reduction  gear  by  a  1 .5-hp 

Boston  Ratiotrol  motor.  The  towing  speed  was  regulated  within  an  accuracy  of 

0.1  percent.  The  system  was  able  to  achieve  towing  speeds  between  1  and 
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300  cm/sec  (Rc  =  2.5x10  -  7.5x10  )  for  the  present  study. 

3 . 2  Wind  Tunnel 

While  the  flow  visualization  experiments  were  conducted  in  the  water  towing 
tank,  a  wind  tunnel  was  used  to  measure  the  lift,  drag  and  pitching  moment  and 
to  survey  the  velocity  field  around  the  delta  wing.  The  venturi  wind  tunnel, 
located  at  the  University  of  Washington's  Department  of  Aeronautics  and 
Astronautics  is  an  open  circuit  type  and  is  6.5  meter  long.  The  test  section 
is  of  a  hexagonal  cross  section,  1  meter  wide,  0.8  meter  high  and  1  meter  long. 
The  wind  tunnel  is  powered  by  a  10  hp  D.C.  motor  driving,  by  direct  coupling,  a 
1.6  meter,  four  bladed  propeller.  The  model  is  mounted  in  the  test  section 
using  a  three-point  support  system  consisting  of  two  fixed  bayonets  forward  and 
a  pitch  arm  aft.  The  maximum  speed  of  the  tunnel  is  32  m/sec  (Rc  =  U^  c/v  = 
8.5x10^)  with  a  clear  test  section.  The  background  turbulence  in  the  tunnel 
is  quite  high  and  exceeds  1  percent  even  with  the  addition  of  several  turbulence 
management  devices  (honeycombs,  screens,  etc.). 

3.3  Delta  Wing  Models 

Four  different  delta  wing  models  were  used  in  the  present  investigation. 

For  the  flow  visualization  experiments  conducted  in  the  water  tank,  a  delta 
wing  with  a  leading  edge  sweep  of  45*  and  a  HACA  0012  profile  at  each  spanwise 
section,  a  delta  wing  with  similar  sweep  angle  but  having  a  sharp  leading  edge 
and  a  flat  surface,  and  a  60*  sweep  wing  with  a  sharp  leading  edge  were  used. 


The  root  chord  of  all  three  wings  was  25  cm.  A  larger  model  having  a  root 
chord  of  40  cm,  a  60*  sweep  angle  and  a  sharp  leading  edge  was  used  in  the 
wind  tunnel  tests.  The  blunt  leading  edge  model  was  made  of  aluminum,  while 
the  sharp  leading  edge  models  were  made  of  Plexiglas. 

3.4  Leading  Edge  Excitation 

To  perturb  the  shar  layer  originating  at  the  leading  edge  of  a  delta  wing, 
a  variety  of  methods  are  available  including  the  following: 

(a)  a  passive  cavity  system  which  is  tuned  to  a  frequency  range  corres¬ 
ponding  to  the  natural  frequency  of  the  vortex  shedding.  The  tuning 
is  provided  by  a  piston  or  a  similar  device  which  changes  the  volume 
of  the  cavity  so  that  its  natural  resonance  frequency  is  comparable  to 
the  shedding  rate. 

(b)  acoustical  means  which  for  example  may  utilize  an  array  of  speakers 
perpendicular  to  the  cross  section  of  the  wing.  The  speakers  are  used 
to  provide  a  velocity  perturbation  along  a  slot  in  the  leading  edge  of 
the  delta  wing. 

(c)  a  piezoelectric  array  embedded  in  the  wing  near  the  leading  edge.  The 
vibration  of  the  piezoelectric  surface  perpendicular  to  the  flow  would 
provide  the  disturbance.  In  a  field  application  this  has  the 
advantage  that  there  is  less  chance  of  fouling  because  no  slots  are 
used . 

(d)  a  heating  element  could  be  utilized  to  provide  an  unsteady  perturbation 
to  the  flow  before  separation.  Because  this  device  may  not  respond  to 
as  high  frequency  as  the  above  methods,  it  may  be  most  useful  in 
hydrofoil  applications  where  the  frequency  requirements  are  lower. 

(e)  blowing  and/or  sucking  through  spanwise  slots  similar  to  those  in 

(b)  above.  This  would  provide  a  larger  amplitude  disturbance  than  the 
acoustical  technique. 

(f)  mechanical  devices,  such  as  a  vibrating  ribbon,  in  the  shape  of  leading 
edge  slots  or  other  configurations. 


KS 


In  this  preliminary  phase  of  the  research,  only  method  (e)  above  was  tested 
Periodic  blowing  or  suction  was  applied  through  both  leading  edges  of  a  delta 
wing  model.  A  slot  0.2-mm  wide  was  formed  at  each  leading  edge  and  plumbed  to 
•  system  of  grooves  inside  the  model.  In  the  towing  tank  experiments,  a  con¬ 
stant  head  tank  was  used  together  with  a  solenoid  value  to  inject  water  at  a 
given  frequency  and  amplitude  out  of  the  leading  edge  slots.  The  solenoid 
valve  was  driven  by  a  square-wave  signal  generator  having  the  desired  frequency 
A  100-psi  compressed  air  line  was  used  to  periodically  inject  air  through  the 
leading  edge  slots  in  the  model  placed  in  the  wind  tunnel.  The  secondary  air 

flow  rate  was  varied  using  a  pressure  regulator  and  the  perturbation  frequency 

* 

was  controlled  using  a  home-made  rotary  valve  driven  by  an  external  motor. 

The  same  valve  was  used  together  with  a  vacuum  chamber  to  achieve  periodic 

suction  in  the  wind  tunnel  experiments. 

The  peak  injection  or  suction  speeds  from  the  leading  edge  slots  varied 

over  the  wide  range  of  0  to  2  U^.  While  the  upper  limit  of  this  speed  may 

seem  excessive,  the  resulting  perturbation  were  still  weak  because  of  the 

extremely  small  area  of  the  leading  edge  slot  compared  to  the  area  of  the  delta 
2  2 

wing  (1/2  p  U.  A.  «  1/2  p  U  A  .  ). 

j  j  «  wing 

3.5  Flow  Visualization  Methods 

The  water  towing  tank  was  used  to  conduct  flow  visualization  experiments  to 
assess  the  nature  of  the  discrete  vortices  shed  from  a  delta  wing  with  and 
without  perturbation  at  its  leading  edge.  Food  color  and  fluorescent  dyes  were 
used  in  the  present  investigation.  The  food  color  dyes  were  illuminated  with 
conventional  flood  lights  and  were  used  to  obtain  overall  views  of  the  flow 
field  using  a  camera  mounted  on  top  of  the  towing  tank.  The  fluorescent  dyes 
were  excited  with  sheets  of  laser  light  projected  in  the  desired  plane.  To 
produce  a  sheet  of  light,  a  5-W  argon-ion  laser  (Spectra  Physics,  Model  164) 
was  used  with  a  mirror  mounted  on  an  optical  scanner  having  a  720-Hz  natural 
frequency  (General  Scanning,  Model  G124).  A  sine-wave  signal  generator,  set  at 
a  frequency  equal  to  the  inverse  of  the  camera  shutter  speed,  drove  the  optical 
•canner  to  produce  light  sheets  approximately  1-mm  thick. 

*The  required  frequencies  in  the  wind  tunnel  experiment  were  higher  than  those 
used  in  the  water  tests. 
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Side  views  of  the  flow  field  were  obtained  using  a  vertical  sheet  of  laser 
light  in  the  x-y  plane  at  z  =  10  cm  (40%  of  the  semispan  of  the  45*  sweep  wing) 
and  a  camera  towed  with  the  wing  but  located  outside  the  tank,  as  shorn  in 
Figure  1.  End  views  were  obtained  using  a  vertical  sheet  of  laser  light  in  the 
y-z  plane  at  x  =  20  cm  (80%  of  the  root  chord)  and  an  underwater  camera  towed 
behind  the  wing  using  the  same  carriage. 

Dye  sheets  or  dye  lines  were  seeped  into  the  boundary  layer  through  a  system 
of  slots  and  holes  on  the  suction  side  of  the  delta  wings.  The  slots  were 
0.2-mm  wide,  milled  at  a  45*  angle  to  minimize  the  flow  disturbance.  The  holes 
were  0.4  mm  in  diameter,  spaced  at  1  cm  center  to  center. 

Dye  was  also  placed  in  the  flow  field  by  laying  several  thin,  horizontal 
sheets  prior  to  towing  the  wing.  The  dye  layers  remained  thin,  about  1-mm 
thick,  due  to  the  inhibition  of  vertical  motion  caused  by  introducing  a  weak 
saline  stratification  in  the  tank.  The  dye  layers  remained  quiescent  until  dis¬ 
turbed  by  the  flow  field  on  and  around  the  wing.  Thus,  the  boundary- layer  flow 
as  well  as  the  potential  flow  could  be  observed  since  the  dye  layers  existed  in 
both  flow  regions.  Cine  films  of  the  runs  presented  in  this  report  are  avail¬ 
able  upon  request. 

-4  4 

For  the  dye-layer  runs,  the  density  gradient  in  the  tank  was  about  10  g/cm  , 
which  yields  a  Brunt-Vaisala  frequency  of  N  =  0.05  Hz.  The  time  scale  asso¬ 
ciated  with  this  weak  stratification  is  large  compared  to  a  typical  convection 
time  scale;  hence,  the  stratification  should  have  little  effect  on  the  dynamics 
of  the  flow.  This  was  verified  experimentally  by  conducting  two  runs  using  the 
dye  injection  method,  in  the  presence  and  in  the  absence  of  stratification  in 
the  tank.  The  visualization  patterns  in  the  two  runs  were  indistinguishable. 

3.6  Velocity  Probes 

Hot-film  and  hot-wire  probes  were  used  to  survey  the  velocity  field  around 
the  delta  wing  mounted  in  the  venturi  wind  tunnel.  Miniature,  single-element 
probes  (TSI  Model  1260)  having  a  diameter  of  0.025  mm  and  a  sensing  length  of 
0.25  mm  were  used.  A  three-dimensional  traverse  was  used  to  survey  the  longitu¬ 
dinal  velocity  fluctuations  in  the  vicinity  of  the  delta  wing.  Conventional 
turbulence  statistics  such  as  the  mean,  root-mean- square,  spectrum,  correlation 
and  probability  distribution  were  computed  from  the  velocity  signals. 
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The  fast-response  probe  signals  were  first  recorded  on  an  instrumentation 
tape  recorder  (Hewlett-Packard  Model  3964A) .  The  FM  analog  tape  recorder  has 
4  channels,  uses  6  mm  tape  and  accepts  very  low  frequencies  including  DC.  The 
data  are  stored  in  the  analog  tape  and  digitized  later  using  the 
analog-to-digital  converter  on  the  NOVA  800  minicomputer.  Data  analysis  are 
then  performed  on  the  NOVA  and  on  the  mainframe  computer  PRIME  P750. 

3.7  Force  Balance 

The  venturi  wind  tunnel  described  in  Section  3.2  has  a  three-component 
balance  system  to  measure  lift,  drag  and  pitching  moment.  The  delta  wing  is 
mounted  on  a  three-point  support  type  system,  consisting  of  two  fixed  bayonets 
forward  and  a  pitch  arm  aft.  The  force  readings  appear  on  a  motor  driven  beam 
balance  having  a  resolution  of  10  gram.  However,  due  to  the  continuous  motion 
of  the  force  balance  caused  by  insufficient  damping  in  the  feedback  circuit, 
the  smallest  force  change  that  was  accurately  readable  was  about  one  order  of 
magnitude  higher  than  the  nominal  resolution.  As  will  be  discussed  in 
Section  8,  the  balance  resolution  was  not  sufficient  to  detect  the  expected 
changes  in  the  aerodynamic  forces  due  to  the  vortex  control  device.  The  short 
period  available  for  the  present  project  (6  months)  did  not  allow  enough  time 
for  considering  alternative  facilities  including  building  a  new  force  balance. 
However,  during  the  second  phase  of  this  research  a  new  force  balance  with 
sufficient  resolution  will  be  designed  and  built. 
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4.  FLOW  VISUALIZATION  RESULTS 

Three  delta  wing  models  were  tested  in  the  18-m  towing  tank.  Different 

views  of  the  flow  field  around  the  perturbed  and  unperturbed  lifting  surfaces 

were  obtained  using  food  color  and  fluorescent  dyes.  Figure  3  shows  a  plane 

view  of  the  60* ,  sharp  leading  edge  delta  wing  at  an  angle  of  attack  of  a  =  10® 

4 

and  a  root-chord  Reynolds  number  of  1.25x10  .  The  flow  is  from  top  to  bottom 

and  the  flow  field  on  the  left  side  of  the  wing  is  visualized  using  a  dye  slot 

at  the  leading  edge.  Note  that  the  dye  sheet  has  alternating  dark  and  light 

regions  near  the  leading  edge  and  extending  along  it.  Although  these  regions 

become  more  obscure  as  they  are  wrapped  up  into  the  large  vortex  structure, 

they  do  indicate  that  a  substructure  is  present  in  the  flow  field. 

Figure  4  is  a  side  view  of  the  45s  sweep  delta  wing  Which  had  a  NACA  0012 

profile  at  each  spanwise  section.  The  angle  of  attack  is  a  =  10®  and  the 

4 

Reynolds  number  is  R  =  2.5x10  .  The  fluorescent  dye  is  illuminated  with  a 
vertical  sheet  of  laser  light  parallel  to  the  flow  direction  at  z  =  10  cm.  In 
this  picture,  the  pairing  process  between  two  discrete  vortices  is  clearly 
indicated. 

A  cross-sectional  view  of  the  discrete  leading  edge  vortices  is  depicted  in 
the  sequence  of  photographs  in  Figure  6.  In  here,  the  60®,  sharp  leading  edge 

4 

delta  wing  at  an  angle  of  attack  of  a  =  10®  and  a  Reynolds  number  of  R  =  1.25x10 

c 

is  visualized  using  fluorescent  dye  seeped  from  the  leading  edge  slot  and 

excited  with  a  vertical  sheet  of  laser  light  perpendicular  to  the  axis  of  the 

primary  leading-edge  vortex.  Six  frames  from  the  cin6  film  are  shown  in 

* 

Figure  6  at  times  t  =  0,  0.14,  0.20,  0.34,  0.49  and  0.56  sec.  It  is  clear 
that  the  dye  rolls  up  into  regions  of  strong  concentration  immediately  after 
separation,  much  the  same  as  in  the  case  of  a  free  shear  layer.  The  discrete 
vortices  are  separated  by  a  very  thin  braid  of  dye  and  two  vortices  are  observed 
to  roll  around  each  other  and  merge  to  form  a  single  larger  vortex.  This  pro¬ 
cess  repeats  itself  at  least  three  times  before  reaching  the  trailing  edge  of 
the  wing  at  the  Reynolds  number  range  used  in  the  present  investigation. 

The  cin6  films  were  used  to  determine  the  natural  shedding  frequency,  fQ, 
by  counting  the  number  of  vortices  passing  a  fixed  position  near  the  leading 
•dge  per  unit  time.  The  shedding  frequency  was  independent  of  the  wing's 

*The  origin  of  time  is  arbitrary. 
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sweep  angle  and  leading  edge  shape.  At  constant  U  ,  the  observed  number  of 

00 

vortices  shed  per  unit  time  tended  to  decrease  as  a  increased.  At  constant  a, 
the  measured  frequency  is  proportional  to  the  square  root  of  the  towing  speed. 
A  least  squares  fit  of  the  data  at  a  =  15*  yields  the  empirical  relation: 

f  =  1.3/u  ;  a  =  15*.  R  =  1 . 25x10*  -  3.33xl05, 

O  v  co  c 


or 


f  c/U 


o  00 


=  1625/ 


If  the  vortices  are  shed  at  a  constant  Strouhal  number,  this  would  imply 
that  the  length  scale  also  increases  as  the  square  root  of  the  velocity. 

Since  the  frequency  is  very  sensitive  to  slight  changes  in  the  shear  layer, 
this  result  may  have  been  influenced  by  increasing  the  dye  flow  rate  at  higher 
velocities  in  order  to  improve  the  visual  image. 

The  principal  result  of  the  visualization  results  depicted  in  Figures  3,  4 
and  6  is  that  the  classical  large  vortices  on  delta  wings  originate  as  a  series 
of  smaller  vortices  shed  from  the  leading  edge  of  the  airfoil.  They  rotate 
around  each  other  and  pair  to  form  larger  vortices  while  simultaneously  moving 
downstream.  As  was  depicted  by  Hoshko  (1976)  for  a  free  shear  layer,  the  pair¬ 
ing  process  is  an  important  contributor  to  the  growth  of  the  vortical  region. 

* 

Turbulent  diffusion  can  correctly  but  imprecisely  describe  that  growth. 

The  best  attack  angle  to  observe  the  shedding  and  pairing  of  the  discrete 
vortices  was  in  the  range  of  10  to  15*.  At  angles  below  5*,  the  pairing  process 
seemed  to  be  inhibited  by  the  proximity  of  the  solid  boundary.  At  large  angles 
of  attack,  the  flow  was  more  turbulent  and  masked  the  observed  phenomenon. 

By  perturbing  the  leading  edge  of  a  delta  wing,  the  shedding  and  the 
pairing  of  the  discrete  vortices  are  modulated  much  the  same  as  the  changes  in 
the  growth  rate  achieved  in  perturbed  free  shear  layers.  The  flow  at  the 
leading  edge  of  the  delta  wing  was  perturbed  in  the  tow  tank  experiments  by 
impulsively  injecting  a  secondary  fluid  from  the  leading  edge  slots.  The 
results  for  a  subharmonic  perturbations  are  shown  in  the  sequence  of  photo¬ 
graphs  in  Figure  7.  The  run  conditions  were  the  same  as  those  depicted  in 
Figure  6,  but  a  square-wave  perturbation  was  applied  with  a  peak  injection 


speed  equal  to  the  towing  speed  (U 


j 


=  5  cm/sec)  and  a  frequency  equal  half 


*0r  molecular  diffusion  in  the  case  of  a  laminar  flow. 


t  =  0  sec 
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the  natural  shedding  frequency  (f  »  0.5  f  =1.1  Hz).  It  is  clear  that  the 
perturbed  shear  layer  is  quite  different  from  the  natural  one.  The  first  pair¬ 
ing  occurs  closer  to  the  leading  edge  and  the  primary  leading  edge  vortex  is 
aore  intense  and  larger  by  about  30%  than  the  one  observed  in  the  unperturbed 

ease. 

Higher  injection  rates  seem  to  be  counter-productive.  As  shown  in  Figure  8, 

when  the  subharmonic  perturbations  were  doubled  in  amplitude  (U,  =  2  U  = 

3  ® 

10  cm/sec),  the  resulting  large  vortex  was  less  organized  and  the  pairing  pro¬ 
cess  was  not  enhanced.  This  was  also  observed,  Figure  9,  when  the  perturbation 

frequency  was  quarter  of  the  natural  shedding  frequency  (f  =  0.25  f  =  0.55  Hz) 

o 

while  the  perturbation  amplitude  remained  high  (U.  =  2  U  =10  cm/sec). 

3  00 

In  summary,  the  visualization  experiments  indicate  that  the  classical  large 
vortices  on  a  delta  wing  originate  as  a  series  of  smaller  vortices  shed  from  the 
leading  edge.  The  pairing  process  between  these  vortices  is  dynamically  signi¬ 
ficant  and  can  be  modulated  by  perturbing  the  leading  edge  of  the  wing.  Weak, 
subharmonic  perturbations  seem  to  enhance  the  pairing  and  result  in  a  more 
organized  large  vortex. 
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Figure  8.  Higher  Injection  Rate  at  the  Leading  Edge  Seems  to  be  Counter 
Productive.  The  Resulting  Large  Vortex  is  Less  Organized. 
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5.  V»T-nrTTV  MKASURPmriS 

The  large  delta  wing  model  was  tested  at  the  venturi  wind  tunnel  to  obtain 
velocity  and  force  measurements.  The  60s  sweep  angle,  sharp  leading  edge  wing 
is  made  of  Plexiglas  and  has  a  root  chord  of  40  cm.  The  root  chord  Reynolds 
number  in  the  experiment  was  varied  in  the  range  of  1x10  to  4x10  .  Miniature, 
single  element  hot-wire  and  hot-film  probes  were  used  to  survey  the  instan¬ 
taneous  longitudinal  velocity  at  different  positions  within  the  leading  edge 
vortices  and  at  different  run  conditions. 

The  coordinate  system  is  shown  in  Figure  1,  where  the  origin  is  at  the  apex 
of  the  delta  wing,  x  is  along  the  root  chord,  y  is  normal  to  the  model  surface 
and  z  is  in  the  spanwise  direction.  Measurements  were  conducted  at  the  attack 
angles  of  10°  and  28°,  at  two  stations,  0.3  c  and  0.6  c,  and  at  different  y’s 
and  z's.  The  complex,  three-dimensional  shear  layer  was  investigated  without 
perturbation,  with  periodic  injection  at  the  leading  edge  and  with  periodic 
suction. 

A  sample  of  the  instantaneous  longitudinal  velocity  signal,  normalized  with 

the  tunnel  speed,  is  shown  in  Figures  10  through  13.  The  unperturbed  shear  layer 

is  surveyed  for  a  Reynolds  number  R  =  1.3xl05  and  angle  of  attack  a  =  28®.  The 

c 

measurements  were  conducted  along  a  line  perpendicular  to  the  x-y  plane  at 
x  =  0.6  c  and  y  =  0.01  c,  where  C  is  the  root  chord.  Four  different  z's  posi¬ 
tion  are  shown  in  the  figures:  z/c  =  0.349;  0.346  (leading  edge);  0.344  and 
0.326.  All  four  positions  are  very  close  to  the  leading  edge,  yet  dramatic 
differences  exist  between  each  profile,  consistent  with  the  presence  of  an 
intense  shear  layer. 

In  Figure  10,  the  negative  spikes  in  the  instantaneous  longitudinal  velocity 

are  consistent  with  the  passing  of  the  discrete  vortices  at  the  fixed  probe 

location.  As  these  vortices  turn,  lower  speed  fluid  is  transported  on  the  high 

speed  side  of  the  shear  layer.  Note  that  the  mean  speed  at  z/c  =  0.349  is 

higher  than  the  tunnel  speed  (U(x,y,z)  >  U  ).  A  well  defined  mean  shedding 

00 

frequency  is  apparent  and  will  be  computed  from  the  spectrum  peak  of  the  raw 
velocity  data. 

The  negative  spikes  are  still  apparent  in  Figure  11  but  are  now  embedded  in 
the  background  turbulence.  The  mean  speed  is  now  lower.  Positive  spikes  are 
observed  as  the  probe  is  traversed  inward  from  the  leading  edge  (Figures  12  and 
13),  and  the  mean  speed  continues  to  drop  indicating  a  strong  spanwise  shear. 
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The  longitudinal  mean  velocity  profiles  at  different  spanwise  locations  is 
shown  in  Figure  14  for  the  same  run  conditions  depicted  in  Figures  10  through 
13.  At  the  centerline  of  the  delta  wing,  z/c  =  0,  the  velocity  near  the  suction 
surface  is  twice  that  of  the  ambient  speed.  It  drops  to  a  low  of  0.6  U  at 

CD 

z/c  «  0.3  just  inward  of  the  leading  edge,  then  rises  quite  fast  and  drops 
again  just  outside  the  leading  edge.  This  strong  shear  region  is  where  the 
discrete  vortices  are  pairing  and  are  dynamically  interacting. 

The  root-mean- square  profiles  obtained  from  the  same  signal  are  shown  in 
Figures  15  and  16.  In  the  first  figure,  J  u2  is  normalized  with  the  tunnel 
speed  U^,  while  in  Figure  16  the  RMS  velocity  fluctuations  are  normalized  with 
the  local  mean  velocity  U.  The  two  peaks  in  the  velocity  fluctuations  correspond 
to  the  positions  of  maximum  shear  (the  zero  crossings  inward  of  the  leading  edge 
in  Figure  14).  The  RMS  is  as  high  as  27%  of  the  local  mean  velocity  at  these 
regions  of  intense  shear. 

Leading  edge  perturbations  could  dramatically  affect  the  intense  shear 
layer.  From  the  visualization  results,  it  was  found  that  subharmonic  perturba¬ 
tions  could  favorably  change  the  growth  of  the  shear  layer  and  enhance  pairing. 

In  Figures  17  and  18,  the  mean  and  rms  profiles  are  depicted  for  three  cases: 

no  perturbation,  subharmonic  injection  (U.  *  1.73  U  )  and  subharmonic  suction 

3  00 

(U.  =  -0.87  U  ),  where  U.  is  the  mean  injection  speed  from  the  leading  edge 
J  00  J  5 

slot.  The  Reynolds  number  in  all  three  cases  is  R  =  1.3x10  ,  the  angle  of 

c 

attack  is  28°,  and  the  spanwise  profiles  are  taken  at  x/c  *  0.3  and  y/c  =  0.003. 
The  leading  edge  perturbations  shift  the  interface  of  the  shear  layer  inward 
towards  the  wing  centerline.  The  maximum  speed  near  the  leading  edge  is  not 
affected,  but  the  minimum  speed  is  dramatically  reduced  by  the  periodic  suction 
or  injection.  Both  the  mean  and  rms  profiles  appear  to  be  more  complex  in  the 
case  of  subharmonic  suction  perturbation.  Mote  that  the  second  peak  in  the  rms 
profiles  is  not  quite  developed  at  this  upstream  station  (x/c  a  0.3  in  Figure  18) 
compared  to  the  further  downstream  location  shown  in  Figure  15  (x/c  *  0.6). 

The  mean  and  rms  profiles  at  four  different  heights  above  the  suction  sur¬ 
face  of  the  delta  wing  are  shown  in  Figure  19  and  20.  The  Reynolds  number  is 
K  »  1.3xl05  and  the  attack  angle  is  a  =  28*.  Subharmonic  perturbations  were 

applied  at  the  leading  edge  of  the  delta  wing  with  U  =  0.73  U  ,  and  the  span- 

j  ® 

wise  profiles  are  taken  at  x/c  *  0.6  and  y/c  =  0.01,  0.03,  0.05  and  0.08.  It 
is  clear  from  these  four  profiles  that  the  shear  layer  is  extremely  complex. 
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Figure  17.  Mean  Velocity  Profiles,  x/c  =  0.3;  R_  =  1.3 
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a  *  28°;  y/c  =  0.003. 
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Simple  statistical  quantities  such  as  the  mean  and  the  root-mean- square  are 
obviously  not  sufficient  to  understand  the  nature  of  this  three-dimensional, 
unsteady  flow  field.  More  complex  statistical  quantities  will  be  presented  in 
the  next  two  sections. 
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6.  PROBABILITY  DISTRIBUTER 

A  useful  statistical  quantity  that  may  shed  some  light  on  the  phenomenon  of 
discrete  vortices  and  the  effects  of  leading  edge  perturbations  on  the  shear 
layer  is  the  probability  density  function  of  the  longitudinal  velocity  fluc¬ 
tuations  (Lumley,  1970),  which  is  a  measure  of  the  relative  amount  of  time 
that  U(t)  spends  at  various  levels.  We  start  with  plots  of  this  function  for 
different  spanwise  locations  near  the  suction  surface  of  the  unperturbed  delta 
wing.  Figure  21  depicts  the  probability  density  function  at  x/c  =  0.6, 
y/c  =  0.005  for  three  spanwise  positions  very  near  the  leading  edge: 
z/c  =  0.349,  0.348  and  0.344.  The  Reynolds  number  is  R  =  1.3x  105  and 
the  angle  of  attack  is  a  =  28°.  In  the  figure,  the  ordinate  is  a  normalized 

histogram  slot.  Thus,  slot  0  is  the  minimum  longitudinal  velocity  recorded  by 
a  hot-wire  probe  at  this  location,  slot  50  is  the  mean  velocity,  and  slot  100 
is  the  maximum  longitudinal  velocity.  The  abscissa  is  a  normalized  probability 
density  function  such  that  the  area  under  the  curve  is  unity.  Just  outside 
the  leading  edge  (Figure  21a  at  z/c  =  0.349),  the  discrete  vortices  manifest 
themselves  as  negative  spikes  in  the  instantaneous  velocity  record  and  the 
probability  density  function  is  biased  towards  the  low  end  of  the  histogram 
slots  (0  to  50).  This  situation  is  reversed  just  inward  of  the  leading  edge 
(z/c  =  0.344  in  Figure  21c),  which  indicates  higher  probability  of  finding 
velocities  higher  than  the  mean  in  the  instantaneous  velocity  records. 

Similar  trends  but  with  a  shifted  spanwise  position  of  the  shear  layer  are 
observed  at  a  higher  normal  position  y/c  =  0.03,  as  shown  in  Figure  22.  In 
here,  four  spanwise  locations  are  recorded:  z/c  =  0.346  (leading  edge  posi¬ 
tion),  0.345,  0.344  and  0.22.  At  this  higher  normal  position,  low  speed 
fluctuations  are  more  probable  even  inward  of  the  leading  edge  (z/c  =  0.344  in 
Figure  22c),  and  this  trend  reverses  further  inward  of  the  leading  edge  as 
shown  in  Figure  22d  at  z/c  =  0.22. 

Perturbing  the  leading  edge  with  subharmonic  injection  of  secondary  fluid 
will  shift  the  shear  layer  and  modulate  the  pairing  process  as  was  indicated 
in  the  mean  and  root-mean-square  velocity  profiles  shown  in  Section  5.  At  a 
fixed  position,  leading  edge  perturbations  should  undoubtly  affect  the 
probability  density  function.  This  is  demonstrated  in  Figures  23a  and  23b  for 

Uj/U^  =  0  and  1.66,  respectively.  The  Reynolds  number  was  R£  =  2.6x10  and 
the  angle  of  attack  was  a  =  10*.  and  the  hot  wire  was  located  at  x/c  =  0.3, 
y/c  =  0.01  and  z/c  =  0.17.  The  probability  density  function  is  wider  in  the 
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Figure  23.  Effects  of  Leading  Edge  Perturbations  on  Probability 

Density  Function,  x/c  «  0.3;  y/c  ■  0.01;  Ec  »  2.6xl05; 
a  «  10*;  z/c  «  0.17 


case  of  subharmonic  leading  edge  perturbations.  The  negative  spikes  are  also 
more  intense  as  indicated  by  the  higher  probability  in  the  lower  speed  regions. 

The  probability  density  function  in  the  case  of  subharmonic  auction  pertur¬ 
bations  is  indicated  at  a  different  run  condition  in  Figure  24.  The  maximum 
auction  speed  from  the  leading  edge  slot  was  U^/U^  =  ~  0.87  and  the 
perturbation  frequency  was  f  =  0.5  f  =30  Hz.  Other  parameters  of  this  run 

5  0 

were  R  =  1.3x10  ,  a  ■  28°,  and  the  probe  position  was  x/c  *  0.3,  y/c  =  0.27 
c 

and  z/c  =  0.17.  At  this  position  the  probability  density  function  is  quite 

narrow  but  is  still  biased  towards  the  low  velocity  fluctuations. 

More  spanwise  variations  of  the  probability  density  function  is  shown  in 

Figures  25,  26  and  27,  at  three  different  heights  above  the  suction  surface  of 

the  delta  wing:  y/c  =  0.01,  0.05  and  0.08,  respectively.  In  all  three  figures 

the  Reynolds  number  was  R  =  1.3x10  ,  the  angle  of  attack  was  a  =  28°, 

c 

and  the  leading  edge  was  perturbed  with  a  subharmonic  perturbation  having  a 

maximum  injection  speed  of  U.  =  1.73  U  .  The  hot  wire  probe  was  always 

J  ® 

at  a  fixed  downstream  distance  of  x  =  0.6  c.  Again,  useful  insight  into  the 
behavior  of  the  leading  edge  vortices  and  the  accompanying  shear  layer  could 
be  derived  from  observing  the  changes  in  the  probability  density  function  with 
height  and  spanwise  location. 
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Figure  27.  Spanwises  Changes  of  Probability  Density  Function. 
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Figure  27.  (Continued) 


7.  SPECTRAL  A»AL.VSTS 


Two  more  statistical  quantities  are  useful  in  analyzing  the  discrete 
vortices  shed  from  the  leading  edge  of  the  delta  wing  and  the  effects  of  the 
perturbation  device  on  these  vortices.  These  are  the  autocorrelation  and 
spectrum  of  the  longitudinal  velocity  fluctuations.  The  autocorrelation 
coefficient  is  the  correlation  between  the  values  of  u  at  two  different  times 
nondimensionalized  using  the  mean  square: 


*(T,  -  Ult)- 

2 

u 


where  the  overbar  denotes  time  average.  This  coefficient  indicates  how  the 
values  of  the  longitudinal  velocity  fluctuations  at  different  times  are 
related.  The  Fourier  transform  of  the  autocorrelation  coefficient  is  the 
power  spectral  density: 
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which  is  the  normalized  spectral  distribution  of  the  squared  intensity  of  the 
longitudinal  velocity  fluctuations  (see,  for  example,  Hinze,  1975;  Tennekes  & 
Lumley,  1972). 

Variations  of  these  two  statistical  quantities  are  shown  in  the  next  four 
figures  for  different  run  conditions.  In  Figure  28,  the  Reynolds  number  was 
R  =  1.3x10  ,  the  angle  of  attack  was  a  =  28°,  and  the  leading  edge  vortices 
were  not  perturbed  (U.  =0).  The  probe  was  located  at  x/c  =  0.6,  y/c  =  0.005 

and  four  different  spanwise  locations:  z/c  =  0.351,  0.349,  0.348  and  0.344. 

The  peak  in  the  velocity  spectrum  and  the  period  of  oscillation  in  the  auto¬ 
correlation  correspond  to  the  natural  shedding  frequency  of  the  leading  edge 
vortices.  At  this  height  very  close  to  the  wing  surface,  that  peak  is  not  very 
strong  but  become  more  pronounced  as  the  probe  moves  inward  from  the  leading 
edge  (Figure  28d).  The  first  zero  crossing  of  the  correlation  coefficient  is  an 
indication  of  the  integral  time  scale.  This  is  about  15  msec  at  z/c  *  0.351 
decreasing  to  about  6  msec  at  z/c  =  0.344. 

The  spectrum  and  autocorrelation  coefficient  at  the  same  run  conditions  but 
at  a  higher  elevation  above  the  wing,  y/c  =  0.03,  are  depicted  in  Figure  29 

for  four  different  spanwise  locations:  z/c  =  0.346,  0.345,  03.44  and  0.22. 
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Figure  29.  Spanwise  Changes  of  Spectrum  and  Autocorrelation 
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At  this  elevation,  the  peak  in  the  spectrum  corresponding  to  the  natural 
shedding  frequency  is  more  pronounced,  particularly  at  z/c  »  0.344  (Pigure  29c). 

The  effects  of  subharmonic  perturbation  at  the  leading  edge  of  the  delta 
wing  are  shown  in  Figures  30  and  31  at  two  different  heights  above  the  wing: 
y/c  >0.03  and  0.08,  respectively.  The  Reynolds  number  and  angle  of  attack  are 
the  same  as  that  in  the  run  depicted  in  Figure  29,  but  a  periodic  injection 
perturbation,  with  maximum  jet  speed  from  the  leading  edge  of  =  1.73  U^, 
is  applied  in  the  runs  shown  in  Figures  30  and  31.  Different  spanwise  loca¬ 
tions  are  depicted  in  the  two  figures.  Comparisons  between  Figures  29b  and  30a, 
and  Figures  29c  and  30b,  which  represent  the  same  position  relative  to  the  wing 
in  the  two  respective  cases  with  and  without  leading  edge  perturbations,  reveal 
differences  in  the  spectrum  and  autocorrelation  that  are  consistent  with  the 
shift  in  the  position  of  the  shear  layer  and  the  corresponding  change  in  the 
mean  velocity  and  the  root-mean-square  (Section  5). 

At  the  higher  elevation  y/c  =  0.08  shown  in  Figure  31,  the  spectrum  peak  is 
weak  near  the  leading  edge,  becoming  more  pronounced  at  z/c  »  0.321  and  0.308 
(Figures  31d  and  31e),  consistent  with  the  general  outline  of  the  shear  layer 
as  sketched  in  Figure  1. 
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Figure  31.  Spanwise  Variations  of  Spectrum  and  Autocorrelation, 
x/c  »  0.6;  y/c/  »  0.08;  Rc  *  1.3xl05;  a  -  28*; 

Uj/U.  «  1.73 


» 


73 


BL53G.  1.  343-363 


CORTBLS36 


CHANNEL  i 


LDG  F 


b.  z/c  -  0.331 


Figure  31.  (Continued) 


74 


C0ATBLS36 


CHANNEL  1 


LOG  F 


J _ 1 _ I _ L 


- 1 - 

30 

SECONDS 


— i — 

40 

10 


— r* 
50 


-3 


eo 


e.  z/c  «  0.308 


Figure  31.  (Continued) 


8.  force  mkasiie  exerts 


The  three-component  force  balance  located  in  the  venturi  wind  tunnel  was 
used  to  measure  lift,  drag  and  pitching  moment  on  the  delta  wing  model.  Three 
force  readings  appear  on  a  motor  driven  beam  balance  having  a  resolution  of 
10  gram,  and  the  readings  are  used  in  a  computer  program  that  takes  into  account 
the  balance  position  with  no  wind  in  the  tunnel  and  the  effects  of  the  model 
support.  The  output  of  the  program  is  the  actual  lift,  drag  and  pitching 
moment  on  the  wing.  For  the  basic  unperturbed  case,  lift,  drag  and  moment 
coefficients  were  determined  for  different  wind  speeds  and  angles  of  attack. 

The  results  were  in  reasonable  agreement  with  the  available  literature,  parti¬ 
cularly  the  exhaustive  set  of  measurements  conducted  by  Wentz  &  Kohlman  (1968) 
for  thirteen  different  delta  wings. 

The  force  and  moment  coefficients  were,  as  expected,  independent  of  Reynolds 
number.  The  sharp  leading  edge  on  the  delta  wing  fixes  the  separation  point 
and  the  vortex  lift  is  mostly  an  inviscid  phenomenon  invarient  to  the  Reynolds 
number,  at  an  attack  angle  of  a  =  10°,  the  three  coefficients  were  C  =0.5, 

L 

C  =  0.1  and  C  =  -0.1.  At  a  =  28°,  the  lift,  drag  and  moment  coefficients  were 
D  M 

C  =  1.15,  C  =  0.6  and  C  =  -0.1. 

L  D  M 

When  the  wing  was  perturbed  at  the  leading  edge  using  subharmonic  suction 
or  injection,  no  measurable  change  of  the  forces  was  observed.  This  result 
constituted  a  setback  to  the  present  project  and  warranted  further  considera¬ 
tion.  The  nominal  resolution  of  the  force  balance  used  in  the  present  experi¬ 
ment  was  10  gram  force.  However,  due  to  the  continuous  motion  of  the  balance 
caused  by  insufficient  damping  in  the  feedback  circuit,  the  needle  moved  con¬ 
tinuously  over  a  range  that  is  about  an  order  of  magnitude  higher  than  the 
nominal  resolution.  Attempts  were  made  to  read  a  mean  position  within  this 
range  of  movement.  Changes  in  the  force  that  are  smaller  than  the  needle 
movement  are  obviously  going  to  be  hard  to  resolve.  At  a  tunnel  speed  of 

U  =  500cm/sec  and  an  attack  angle  of  a  =  10°  the  lift  force  was  approximately 
00 

70  gram.  If  the  leading  edge  perturbation  device  that  is  being  developed  in 

* 

the  present  project  is  to  change  the  lift  coefficient  by  10%,  the  change  in 
the  balance  reading  would  be  only  7  gram;  obviously  too  small  to  resolve  with 
the  present  force  measuring  device.  Even  at  the  higher 


*This  is  probably  a  significant  change  for  an  actual  fighter  aircraft. 


speed  of  U  *  1000  cm/sec  and  the  higher  angle  of  attack  of  a  =  28s,  the 
00 

resulting  lift  force  is  slightly  higher  than  600  gram.  A  10%  change  translates 
to  about  60  gram  change  in  the  reading  of  the  balance.  While  this  is  higher 
than  the  nominal  resolution  (10  gram),  this  reading  is  within  the  needle  fluc¬ 
tuations  due  to  the  insufficient  damping  in  the  feedback  circuit  as  noted 
above . 

The  short  period  available  for  the  present  project  (6  months)  did  not 
allow  enough  time  for  considering  an  alternative  force  balance,  particularly 
building  a  new  one.  However,  during  the  second  phase  of  this  research  a  new 
force  balance  with  sufficient  resolution  will  be  designed  and  built. 


9.  SUMMARY 

Based  on  our  recent  finding  that  a  leading  edge  vortex  on  a  delta  wing  in 
steady  flight  consists  of  a  series  of  discrete  smaller  vortices,  a  device  was 
proposed  to  modulate  these  discrete  vortices  to  favorably  affect  the  flow 
around  the  wing.  The  purpose  of  the  present  phase  of  the  research  reported 
herein  is  to  establish  the  proof  of  concept  and  to  understand  the  complex  flow 
field  under  consideration.  To  this  end,  the  perturbed  and  unperturbed  flow 
fields  around  a  delta  wing  at  constant  angle  of  attack  were  studied  experi¬ 
mentally  in  a  water  towing  tank  and  a  high  speed  wind  tunnel.  Flow  visualiza¬ 
tion,  fast-response  velocity  probe  surveys,  as  well  as  force  measurements  were 
conducted  to  study  the  shear  layer  formed  on  the  suction  side  of  the  wing,  and 
to  determine  the  effects  of  periodic  suction  or  injection  from  a  leading  edge 
slot  on  the  shedding  and  pairing  of  the  discrete  vortices. 

The  flow  visualization  experiments  utilized  food  color  and  fluorescent 
dyes  illuminated  with  flood  lights  or  sheets  of  laser  projected  in  the  desired 
plane.  Periodic  injection  from  a  leading  edge  slot  significantly  changed  the 
shear  layer  on  the  suction  side  of  the  wing.  Subharmonic  perturbation  with 
relatively  low  amplitude  seemed  to  enhance  the  pairing  process  and  result  in  a 
more  organized  large  vortex. 

Miniature  hot-wire  and  hot-film  probes  were  used  to  survey  the  instan¬ 
taneous  longituding  velocity  at  different  positions  within  the  leading  edge 
vortices.  Different  statistical  quantities  were  computed  including  the  mean, 
root-mean- square,  probability  density  function,  autocorrelation  and  spectrum. 
The  instantaneous  signals  clearly  showed  positive  and  negative  velocity  spikes 
consistent  with  the  passing  of  the  shed  vortices  by  the  fixed  probe.  Leading 
edge  perturbations  dramatically  affected  the  intense  shear  layer  as  evidenced 
by  the  marked  changes  in  the  different  statistical  quantities. 

The  expected  changes  in  the  aerodynamic  forces  were  not  observed  in  this 
preliminary  phase  of  the  research.  Careful  analysis  of  the  force  balance  used 
in  the  present  experiments  indicated  that  these  changes  are  most  probably 
within  the  accuracy  of  the  balance  system. 

The  results  of  Phase  I  research  indicate  that  the  discrete  vortices  could 
be  strongly  modulated  with  periodic  suction  or  injection  from  a  leading  edge 
slot.  In  the  next  phase,  a  force  balance  with  adequate  resolution  will  be 
designed  and  built,  and  the  perturbation  device  will  be  optimized  for  field 
applications . 
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